Stripe Charge Ordering in Triangular-Lattice Systems 
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We investigate the ground-state properties of a t2g-orbital Hubbard model on a triangular lattice 
at electron density 5.5 by using numerical techniques. There appear several types of paramagnetic 
phases, but we observe in common that one or two orbitals among three orbitals become relevant 
due to the effect of orbital arrangement. It is found that charge stripes stabilized by the nearest- 
neighbor Coulomb interaction consist of antiferromagnetic/ferro-orbital chains for small Hund's 
couphng, while there occurs stripe charge ordering with ferromagnetic/antiferro-orbital chains for 
large Hund's coupling. 
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In recent years, a wide variety of exotic phenomena 
in layered cobalt oxides have attracted much attention 
in the research field of condensed-matter physics. In 
particular, the discovery of superconductivity in water- 
intercalated cobaltite Nao.35CoO2-l.3H2O [l| has trig- 
gered current high activity on the study of supercon- 
ductivity in frustrated triangular-lattice systems. On the 
other hand, magnetic transitions have also been observed 
in Naa;Co02 [Midi' In such materials, the competition be- 
tween various types of ordered states should arise due to 
geometrical frustration. In addition, it has been stressed 
that the i2g-orbital degree of freedom plays a crucial 
role in causing the rich phase diagram in the cobaltites. 
In fact, multiorbital models have been investigated to 
understand the mechanism of exotic superconductivity 
and novel magnetism in triangular-lattice cobalt oxides 
0, H, 0, Q]. To gain a deep insight into the behavior 
of such complex cobalt-based systems, it is important 
to clarify the characteristics of spin-charge-orbital states 
taking account of strong electron correlations 

In this paper, we study the ground-state properties of 
the t2g-orbital Hubbard model on a triangular lattice at 
electron density 5.5. The Hamiltonian is given by 
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FIG. 1: (a) Triangular lattice with three directions depicted 
by arrows. An eight-site cluster enclosed by dotted lines is 
considered in our numerical calculations, (b) Schematic view 
of level splitting among t2g orbitals. 



Fig. 1). The hopping amplitudes are given by t^y^^y = 
[ddn), tyz^yz ~ ^yz,zx ~ ^zx,yz ~ ^zx.zx ~ {ddn) 
for the X direction, t^y = {9/16){dda) + (l/4)(c?d7r), 
f^z^yz = (3/4)(dd7r), i^,,,, - (l/4)(dd7r), t^,,,, = 
^zx,yz = {V^/4)iddn) for the u direction, and t^ ^ = t^_^, 
except for ty^ .,^ = t^^^yz = (— •\/3/4)((id7r) for the v di- 
rection. Note that there is no hopping between the aig 
(xy) and e'„ {yz and zx) orbitals due to their different 
symmetry. Hereafter, we take {dd'K)=l as the energy 
unit. A is the level splitting between the oig and Cg 
orbitals, and V is the nearest-neighbor Coulomb inter- 
action. In the local interactions, U, U', J, and J' are 
intraorbital Coulomb, interorbital Coulomb, exchange, 
and pair hopping interactions, respectively. Note the 
well-known relations J=J' and U=U'+J+J' . We ana- 
lyze the model for eight-site {N=8) systems with periodic 
boundary conditions by the Lanczos diagonalization. In 
this paper, we set {dda)=2, A=6, and U'=20, and inves- 
tigate the dependence on J and V. The dependence on 
{dda), A, and U' will be discussed elsewhere in future. 

Our result is summarized in Fig. 2(a), which is the 
ground-state phase diagram in the {V, J) plane. Com- 
paring the lowest energies in the subspaces of the z com- 
ponent of the total spin, we determine the total spin in 
the ground state. Then, we find that for a large set of pa- 
rameters, the ground state is a paramagnetic (PM) phase 
characterized by total spin zero, but there appear several 
types of PM phases in terms of spin-charge-orbital con- 
figurations. Among them, the phase I is characterized 
by a stripe, in which ferromagnetic (FM) spin chains are 
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FIG. 2: (a) Ground-state phase diagram in the {V, J) plane 
for A=6 and U'—20 at electron density 5.5. Inset denotes 
schematic view of the charge stripe and spin configuration 
in each phase. The size of circle indicates the degree of the 
charge density, (b) Electron densities in the xy, yz, and zx 
orbitals. (c) Electron density in the a orbital, represented by 
the linear combination of the yz and zx orbitals. 

coupled with each other in an antiferromagnetic (AFM) 
manner, as shown in the inset of phase I. Note that 
the charge is uniformly distributed in the background 
of ferro-orbital (FO) arrangement. 

On the other hand, the phases II and III exhibit charge 
order. To clarify what types of charge configurations 
emerge in the present system, we investigate the charge 
structure factor. When we examine the V dependence 
of the charge structure factor for J=0, it is clearly ob- 
served that it develops abruptly around at V—IQ, in- 
dicating the stabilization of a specific charge structure. 
The results suggest a charge stripe composed of one- 
dimensional chains, as shown in the inset of phase II. 
Namely, the spin frustration is removed due to charge 
ordering. It should be noted that a similar stripe charge 
pattern has been suggested for a C0O2 triangular-lattice 
model by Hartree-Fock calculations |3 ■ Note that in the 
phase II, two of three orbitals are found to be fully occu- 
pied due to the effect of FO arrangement, similar to the 
phase I. 

Now we consider the effect of J on the spin structure 
in the stripe charge state. For this purpose, it is useful to 
measure the spin structure factor. For V=12, we observe 
a transition in the spin structure factor at around J=0.6, 
but the stripe charge configuration is not changed. For 
J < 0.6, we find AFM correlations in each one-dimensional 
charge stripe, as shown in the inset of phase II. On the 
other hand, for J>0.6, the spin configuration changes 



to a FM one, as shown in the inset of phase III. Weak 
antiferro-orbital (AFO) correlations exist in the chains, 
since the electrons can gain kinetic energy due to the 
so-called double-exchange mechanism. 

In order to identify the orbital arrangement, it is conve- 
nient to evaluate the electron densities, = '^{{pi-y) / N , 
in the xy, yz, and zx orbitals. In Fig. 2(b), we show n-y 
vs. J for V=12. It is observed that the lower xy orbital is 
fully occupied irrespective of J, while the electron densi- 
ties in the yz and zx orbitals change suddenly at around 
J=0.6. For J<0.6, electrons favorably occupy the zx or- 
bital, but the zx orbital is not fully occupied. On the 
other hand, the yz and zx orbitals are equally occupied 
for J>0.6. 

However, in order to determine the orbital structure, 
we need to pay due attention to the definition of orbitals. 
Here, by using an angle 9i to characterize the orbital 
shape, we introduce operators for a and b orbitals such as 
diaa = e''^'/'^[cos{9i/2)di,y^_,c, + sin{9i/2)dizx,a] and diba = 
e'''/'[-sin(ei/2)di,j,,,, + cos(0i/2)di,,,,<,]. In Fig. 2(c), 
we show the electron density in the a orbital with 9i=6. 
We find that for J=0, the a orbital of 6'=47r/3 is fully 
occupied. Namely, only one orbital becomes relevant for 
small J. This FO arrangement is also found in the phase 
I. On the other hand, for J=l, the electron density does 
not change even when 9 is varied, indicating that two 
orbitals remain active. 

Note that we have suggested the stripe charge order 
from the 8-site calculations, but the shape of the stripe 
may change from straight to zigzag. To confirm this 
point, it is necessary to perform the calculations in larger- 
sized systems, for instance, using a model including only 
e'o- orbitals. This is among possible future developments. 

In summary, in the FM phase, one or two orbitals 
among the t2g orbitals become relevant due to the ef- 
fect of orbital arrangement. The charge stripe stabilized 
by V consists of AFM/FO chains for small J, while it 
changes to FM/AFO ones with increasing J. 
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